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Abstract 

Backgrounds: Tension-type headache is a common psychosomatic disease. However, diurnal variation of headache 
is yet to be clarified, perhaps due to the lack of an appropriate method to investigate it. Like other painful diseases, 
it would be helpful to know if there is diurnal variation in tension-type headaches, both for managing headaches 
and understanding their pathophysiology. The aim of this study was to determine if there is diurnal variation in the 
intensity and exacerbation of tension-type headache. 

Methods: Patients (N = 31) with tension-type headache recorded for one week their momentary headache intensity 
several times a day and their acute headache exacerbations using a watch-type computer as an electronic diary 
(computerized ecological momentary assessment). Multilevel modeling was used to test the effects of time of day 
on momentary headache intensity and on the occurrence of acute exacerbations. 

Results: A significant diurnal variation in momentary headache intensity was shown (P = 0.0005), with the weakest 
headaches in the morning and a peak in the late afternoon. A between-individual difference in the diurnal pattern 
was suggested. On-demand medication use was associated with a different diurnal pattern (P = 0.025), suggesting 
that headache intensity decreases earlier in the evening in subjects who used on-demand medication, while 
headache subtype, prophylactic medication use, and sex were not associated with the difference. The occurrence 
of acute headache exacerbation also showed a significant diurnal variation, with a peak after noon (P = 0.001 5). 

Conclusions: Tension-type headache was shown to have a significant diurnal variation. The relation to 
pathophysiology and psychosocial aspects needs to be further explored. 
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Backgrounds 

Tension-type headache (TTH) is one of the common 
and major diseases in psychosomatic medicine because 
of its high prevalence and supposed association with 
psychosocial factors. Despite its familiarity, diurnal vari- 
ation of headache, which is thought to be a quite basic 
characteristic, is still unclear. One of the reasons was lack 
of an appropriate method to investigate a detailed diurnal 
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variation. One cross-sectional study reported that almost 
half of the TTH patients had their usual onset of head- 
ache attacks during the day [1]. However, this report re- 
lied on the patients' recalled perception regarding what 
time their headaches most frequently had occurred. It 
has often been pointed out that recall bias is inevitable in 
surveys based on subjects' recall. No prospective study 
has been conducted to confirm temporal distribution of 
TTH attacks throughout the day. Moreover, diurnal vari- 
ation of headache intensity in TTH has rarely been inves- 
tigated, especially by prospective studies. 

Recently, ecological momentary assessment (EMA), 
which avoids recall bias, has been used in the area of psy- 
chosomatic medicine to document subjective symptoms 
in subjects' daily lives [2]. EMA is a method to assess and 
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record symptoms at the moment they happen via record- 
ing throughout their daily lives. It is done using paper- 
and-pencil diaries or electronic diaries (i.e. computerized 
EMA). In TTH, consistency between recalled headache 
intensity and momentary headache intensity recorded by 
computerized EMA was shown to be low in our previous 
study [3]. In addition, compared with EMA using paper- 
and-pencil diaries, computerized EMA has the advantage 
of avoiding 'faked compliance', i.e. disguised compliance 
by recording data at times other than those designated, 
by recording the time of data input [4]. This feature is 
particularly appropriate for investigating temporal symp- 
tom patterns. 

Diurnal variation of pain intensity has been discussed 
in several types of pain such as cancer and arthritis pain 
in the context of optimization of pain control and the 
relation to pathophysiology [5,6]. In TTH as well, inves- 
tigation of diurnal variation of headache intensity would 
be helpful for patient education, headache management, 
and the understanding of their pathophysiology. There- 
fore, the aims of this study were to prospectively investi- 
gate the diurnal pattern of TTH in the daily setting 
using computerized EMA and to determine if there is 
diurnal variation of headache intensity in TTH and if 
there is a specific time of day when acute headache 
exacerbations are more likely to occur. 

Methods 

All procedures and materials were approved by the eth- 
ics committee of the University of Tokyo. 

Subjects 

The subjects were participants in a trial of relaxation 
therapy for TTH and overlapped with the subjects in a 
previously reported study [7]. Recruitment was by an ad- 
vertisement on websites from March 2003 to August 
2004. Patients who applied were interviewed and 
screened by well-trained physicians. 

Inclusion criteria included age 20 years or older but 
less than 60 years with any type of TTH according to 
the criteria of the International Headache Society [8] 
and with at least one headache episode per week, on 
average. The patients were excluded if they had head- 
ache other than TTH according to the criteria of the 
International Headache Society [8]; current psychiatric 
disease; history of paranoia or schizophrenia; severe 
physical illnesses; if they had participated in relaxation 
therapy; and if they were employed as a shift worker. 
Diagnosis of psychiatric disease was made according to 
the criteria in the Diagnostic and Statistical Manual of 
Mental Disorders, Fourth Edition, Text Revision [9]. 

Eighty-four subjects applied to participate in the study. 
Ultimately, 54 of 59 eligible subjects were enrolled in the 
relaxation therapy trial; Five declined to participate due 



to scheduling conflicts. After the exclusion of 19 subjects 
with migraine and one subject who worked as a shift 
worker, 34 remained eligible for study. All the subjects 
gave their written informed consent to participate. 

Measurement of momentary headache intensity 

Momentary headache intensity was recorded with 
watch-type computers (Ruputer ECOLOG; 42 grams, 
Seiko Instruments Inc., Tokyo, Japan) used as electronic 
diaries [3,7]. The computer had a screen measuring 
20 x 30 mm and a joystick and button as input devices. 
The subjects were fully instructed how to use the device 
and given manuals before the beginning of the study 
period. They also practiced manipulating the device with 
one of the authors (HK) until they became accustomed 
to its use. 

The subjects wore the watch-type computers for seven 
consecutive days before relaxation therapy started. 
Signal-contingent recordings, which were defined as 
recordings prompted by a beep [2], were programmed to 
occur randomly within an interval of 36 minutes at ap- 
proximately 6:00, 12:00, 18:00, and 24:00. The subjects 
were allowed to postpone input for 30 minutes. Record- 
ings not made within 30 minutes were cancelled. The 
subjects were also asked to record their headache inten- 
sities when waking up and going to bed by choosing 
'waking up' or 'going to bed' from the menu. After 
selecting 'going to bed) computers suspended the beeps 
for signal-contingent recordings until a 'waking up' 
recording was selected to avoid disturbing the subject's 
sleep. Signal-contingent recordings and recordings when 
waking up and going to bed were treated as scheduled 
recordings. 

Event-contingent recordings were defined as record- 
ings that were initiated by the subjects themselves when 
a particular event occurred [2]. In this study, the sub- 
jects were asked to make a recording as soon as possible 
every time their headache became exacerbated as an 
event-contingent recording. 

In both scheduled and event-contingent recordings, 
momentary headache intensity was rated according to a 
visual analog scale (VAS) from 0 to 100 displayed on the 
screen. The words "headache intensity" were displayed 
with the VAS as a question, and the anchor words 
"none" and "most intense" were displayed at the respect- 
ive ends of the scale. The subjects manipulated the joy- 
stick and adjusted the length of the bar so that it 
corresponded to their headache intensity at that 
moment. 

Statistical analysis 

Multilevel modeling was used to investigate diurnal 
variation of headache intensity and acute headache 
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exacerbations, using SAS Proc Mixed and Proc Nlmixed 
(SAS 9.1, SAS Institute Inc., Cary, NC). 

First, time of day was divided into three-hour blocks 
(e.g. 0:00-3:00, 3:00-6:00, etc.) and treated as a categor- 
ical variable. Momentary headache intensity was treated 
as the dependent variable and time of day as the pre- 
dictor. The effect of time was modeled either as fixed or 
random. Therefore, the following models were tested: 
model Al, an unconditional model (no predictor); model 
A2, a model with the effect of time as a fixed effect; 
model A3, a model with the effect of time as a random 
effect. The level 1 intercept was the individual true value 
of the momentary headache intensity when the time of 
day was in the reference block (21:00-24:00), and it was 
modeled as a random effect. The variance-covariance 
matrix (G matrix) was modeled as variance component 
structure. 

Second, the effect of time on the logarithm of the aver- 
age number of headache exacerbations in three hours 
was tested by using a multilevel Poisson model in order 
to investigate what time of day headache exacerbations 
were most likely to occur. The number of headache 
exacerbations for every three hours were calculated for 
each patient as the outcome variables, and time of day 
was divided into three-hour blocks as described above as 
a categorical variable. The time block 3:00-6:00 was 
excluded from the analysis because no subjects had 
headache exacerbation during that time block. Effects 
of time were modeled as fixed. The following models 
were tested: model Bl, an unconditional model (no pre- 
dictor); model B2, a model with the effect of time 
(three-hour blocks) as a fixed effect. The level 1 inter- 
cept (individual true value of the logarithm of the aver- 
age number of headache exacerbations when the time 
of day was in the reference block (21:00-24:00)) was 
modeled as a random effect. 

Third, we investigated the difference of the diurnal 
pattern between chronic and episodic tension-type head- 
ache by analyzing the models with a variable of headache 
subtype and the interaction term of time and headache 
subtype in order to investigate the heterogeneity of diur- 
nal pattern. Headache diagnosed as "headache of 
tension-type not fulfilling the criteria of episodic or 
chronic tension-type headache" was coded as either 
chronic or episodic according to average headache fre- 
quency (> 15 days/month or not). We also investigated 
the effect of medication use on the diurnal variation of 
TTH by analyzing the models with a dichotomous vari- 
able of medication use and the interaction term of time 
and medication use. We tested two types of medication 
use: The use of prophylactic medication and the use of 
on-demand analgesic medication during the study 
period. In addition, we investigated the difference of the 
diurnal pattern between women and men. 



Goodness of fit was compared using a -2 log likelihood 
function and x test when one model was nested in the 
other; otherwise Akaike's Information Criterion (AIC) 
was used. Post-hoc multiple comparisons were con- 
ducted with either Turkey-Kramer correction or Bonfer- 
roni correction. The significance level was set at 0.05. 

Results 

Patient characteristics 

Thirty-four subjects were enrolled in the study and three 
were excluded from further analysis because they were 
unable to complete their recordings for seven days due 
to problems with the computers. Finally, data from 31 
subjects (nine men and 22 women) were analyzed. The 
profiles of all subjects are shown in Table 1, with women 
and men listed separately. The mean age of the subjects 
was 38.4 years (SD, 10.4 years; range, 25-60 years). Five 
subjects had episodic TTH, 24 had chronic TTH, and 
two had headache of tension-type not fulfilling the cri- 
teria of episodic or chronic tension-type headache. 

Table 1 Demographic and medical characteristics of the 
subjects 



Age 

Mean in years (SD) 
Subtype of tension-type headache 
Episodic 
Chronic 
Other* 

Disease duration 
Median in years (range) 

Prophylactic medication 
With 

anxiolytics 

muscle relaxants 

Chinese herbal medicine 

antidepressants 
Without 

Use of on-demand medication 
With 

analgesics (either as a 
single or compound 
preparation) 

anxiolytics 

muscle relaxants 

Without 



Total Women Men 
(N = 31) (N = 22) (N = 9) 



38.4 (10.4) 38.1(11.1) 38.7(9.0) 

5 (16%) 4 1 

24 (77%) 16 8 

2 (6%) 2 0 

7.0 7.5 3.0 

(0.1-24.0) (0.1-24.0) (1.3-23.0) 

13 (42%) 8 5 

9 (29%) 5 4 

5 (16%) 2 3 

4 (1 3%) 4 0 
3(10%) 1 2 

18(58%) 14 4 

19(61%) 15 4 

14(45%) 10 4 

5 (16%) 5 0 
2 (6%) 2 0 

12(39%) 7 5 



SD, standard deviation. 

* headache of tension-type not fulfilling the criteria of episodic or chronic 
tension-type headache. 
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Thirteen subjects answered that they regularly took 
prophylactic medication and 19 answered that they usu- 
ally used analgesics on-demand. The subjects who were 
excluded from the analysis were all men, and two of 
whom had chronic TTH and the other had TTH not ful- 
filling the criteria of episodic or chronic tension-type 
headache. Neither their ages nor headache intensities 
rated by weekly recall were significantly different from 
those of the subjects who completed their recordings. 

Recording profiles 

For all subjects, there were 998 scheduled recordings, 
consisting of 570 signal-contingent recordings, 213 
recordings on awakening, and 215 recordings at bed- 
time. The mean compliance rate for signal-contingent 
recording was 96%. The mean number of scheduled 
recordings was 32.6 per subject. Twenty-three subjects 
made 105 event-contingent recordings and the other 
eight subjects made no event-contingent recordings. The 
number of recordings made in each three-hour time 
block and the number of subjects who made recordings 
in each three-hour time block are shown in Figure 1. 

Diurnal variation of momentary headache intensity 

Of the three models, the model with time as a random 
effect (model A3) had the best fit. The effect of time of 
day was significant in model A3 (F(7, 193) = 3.91, 
P = 0.0005). Multiple comparisons of all pairs of esti- 
mated momentary headache intensity from each time 
block (with Turkey-Kramer correction) revealed that 
momentary headache intensity from 3:00-6:00 was sig- 
nificantly lower than from 15:00-18:00 (t(193) = -3.38, 
adjusted P= 0.019) and 18:00-21:00 (t(193) = -3.35, 
adjusted P = 0.021), and that momentary headache inten- 
sity from 6:00-9:00 was significantly lower than from 



15:00-18:00 (t(193) = -3.73, adjusted P = 0.006) and 
18:00-21:00 (t(193) = -3.78, adjusted P = 0.0051). Model 
A3 is expressed as follows: 
Level 1 equation: 



Yy = n 0i + TZiiTIMEOy + 7r 2i TiME3, y 
+ 7r 7 ,TIME18,7 + cu 



Level 2 equations: 



Tim 



Yoo 
Yko 



Co, 

C kl (k : 



1,2, 



where Yy is each momentary headache intensity for the 

i'th patient; TIMEOy, TIME3y, , and TIME18y are 

dummy variables indicating that the time of day was in 

0:00-3:00, 3:00-6:00, ,or 18:00-21:00, respectively (e.g. 

TIMEOy = 1 if the time of day was in 0:00-3:00, otherwise 
TIME0y = 0). n 0l - is the individual i's true value of moment- 
ary headache intensity when time of day was in the refer- 
ence block (21:00-24:00) and n u , and n 7l are 

individual i's true differences in momentary headache inten- 
sity between corresponding time blocks and the reference 
time block, ys are group means of corresponding ns. Eg, (, oi 
and ( ki (k= 1, 2, , 7) are residuals at each level and in- 
cluding them in the equations means that the intercept and 
the effect of time are modeled as random. The estimated 
momentary headache intensity is plotted in Figure 1. 

G matrix is modeled as a variance component struc- 
ture in this model, which means that ij^s (k = 1, 2, , 

7) are assumed to have the same value of variance. The 
estimated variance of (i <; s was 38.5, which is significantly 
different from zero (P< 0.0001). 



45 



m 35 

! 

a 30 



S 25 



20 




0 3 6 9 12 15 18 21 24 

Number of recordings/ 172/27 41/19 171/30 151/31 149/31 122/31 156/31 141/31 
Number of subjects 

who made recording(s) Time of Day 

Figure 1 Diurnal variation in headache intensity. The estimated mean headache intensity is shown with its standard error for each three-hour 
block, a P<0.05 vs. 3:00-6:00; b P<0.05 vs. 6:00-9:00, using Tukey's multiple comparison. 



Kikuchi ef al. BioPsychoSocial Medicine 201 2, 6:1 8 
http://www.bpsmedicine.eom/content/6/1/1 8 



Page 5 of 7 



Average number of headache exacerbations throughout 
the day 

The model with the effect of time (model B2) fit better 
than the unconditional model, and the effect of time was 
significant (F(6, 30) = 4.81, P = 0.0015). Model B2 is 
expressed as follows: 
Level 1 equation: 



Log(Ave,y) 



n 0i + nuTIMEOij 

7T 6i TIME18,y 



772iTIME6,y 



Level 2 equations: 



Too + (ot 
y k0 (k= 1,2, 



,6) 



The number of headache exacerbations in a three- 
hour block ~ Poisson (Ave,-,) 

where Ave,, is the parameter of the Poisson distribution 
(average number of headache exacerbations per three 
hours) for each time block for the ith patient; and 
TIMEOy, TIME6y, , and TIME18, y are dummy vari- 
ables for the time of day used in the same way as they 
are used in model A3. tt 0 ; is the individual i's logarithm 
of the average number of headache exacerbations when 
time of day was in the reference block (21:00-24:00) and 
tti„ , and tt 6i are individual i's differences in the loga- 
rithmic average number of headache exacerbations be- 
tween corresponding time blocks, and the reference time 
block, ys are the group mean of corresponding us. (o; is 
the residual. The estimated average number of headache 
exacerbations per three hours is plotted in Figure 2. 

Multiple comparisons of all pairs of estimated loga- 
rithmic average number of headache exacerbations 
for each time block (alpha level set at 0.0033 with 



Bonferroni correction) were conducted, and it was 
revealed that the logarithm of the average number of 
headache exacerbations from 0:00-3:00 was significandy 
smaller than from 12:00-15:00 (t(30) = -3.59, P= 0.0012) 
and 15:00-18:00 (t(30) = -3.41, P = 0.0019). 

Effect of headache subtypes on diurnal variation of 
momentary headache intensity and on the average 
number of headache exacerbations throughout the day 

The models with the interaction of time and headache 
subtype did not fit better than the models without it. 
Furthermore the interaction was not significant in the 
models with it. 

Effect of medication use on diurnal variation of 
momentary headache intensity and on the average 
number of headache exacerbations throughout the day 

The models with the interaction of time and prophylac- 
tic medication use did not fit better than those without 
it, either in the analyses of diurnal variation in moment- 
ary headache intensity or in the analysis of diurnal vari- 
ation of headache attacks. Furthermore, the interaction 
of time and prophylactic medication use was not signifi- 
cant in the models including them. 

The model with the effect of time as a random effect and 
the interaction of time and on-demand analgesic medica- 
tion use fit better than the model without the interaction, 
and the interaction was significant (F(7,186) = 2.36, 
P = 0.025) in the analysis of diurnal variation of momentary 
headache intensity. The estimated momentary headache in- 
tensity for each group (patients with and without on- 
demand analgesic medication use) is plotted in Figure 3. 
The model with the interaction of time and on-demand an- 
algesic medication use did not fit better than the model 
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Figure 2 Probability of headache exacerbation. The estimated average number of headache exacerbations per three hours is shown with its 
standard error for each three-hour block. 
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Figure 3 Diurnal variation in headache intensity and on-demand analgesic medication use. The estimated mean headache intensity is 
shown as a solid line for the patients who used on-demand analgesic medication during the study period and as a dashed line for those who 
did not, both with the standard error for each three-hour block. 



without them in the analysis of diurnal variation of head- 
ache attacks. Furthermore the interactions were not signifi- 
cant in the models with them. 

Effect of sex on diurnal variation of momentary headache 
intensity and on the average number of headache 
exacerbations throughout the day 

The models with the interaction of time and sex did not 
fit better than those without it, either in the analyses of 
diurnal variation in momentary headache intensity or in 
the analysis of diurnal variation of headache attacks. Fur- 
thermore, the interaction of time and sex was not sig- 
nificant in the models including it. 

Discussion 

A significant diurnal variation in headache intensity was 
shown in TTH through utilization of computerized 
EMA, with headache being weakest in the morning, 
worsening toward the evening, and peaking afterward. 

When the effect of time was modeled as a random ef- 
fect, the models fit better than when the effect was mod- 
eled as a fixed effect. In addition, the estimated variances 
of those random effects were significantly different from 
zero. These results suggested that there was significant 
between-individual variation in the effect of time. In 
other words, the pattern of diurnal variation of TTH was 
different among patients, while a certain diurnal pattern 
was still observed when averaged for all patients. 

The heterogeneity of diurnal patterns of headache inten- 
sity may stem from different subtypes of tension-type 
headache and sex. However, the analysis failed to show a 
significant difference in diurnal patterns between chronic 
and episodic tension-type headache or between women 
and men. Another possibility is that medication use 



influenced the time course of headache. While prophylac- 
tic medication use did not show any significant difference 
in the diurnal pattern of headache intensity, on-demand 
medication use was related to a significant difference in 
diurnal pattern of headache intensity. The headache inten- 
sity decreased earlier in the evening in patients who used 
on-demand medication during the study period, which 
might reflect the effect of the medication. 

The circadian variation in the probability of acute head- 
ache exacerbations also showed a specific diurnal pattern. 
The estimated average number of headache exacerbations 
was small in the morning, reached a peak just after noon, 
and then gradually decreased towards night. These results 
were consistent with a previous study in which almost half 
of TTH patients reported that their usual headache attack 
onset was during the day [1]. In this prospective study 
using computerized EMA, it was possible to determine 
more accurate and detailed information on the temporal 
distribution of headache exacerbations. Neither subtype of 
headache nor medication use was related to a significant 
difference of circadian variation on the probability of acute 
headache exacerbation. 

Circadian pain intensity variation generally has been 
discussed in relation to endogenous opioid peptides and 
melatonin [5,6]. However, their involvement in TTH has 
not yet been established [10,11]. Furthermore, a previous 
study reported that the diurnal pattern itself was differ- 
ent for plasma melatonin between patients with chronic 
TTH and healthy controls [12]. Further investigation is 
necessary in order to clarify the biological origin of the 
circadian variation of pain intensity in TTH. In addition, 
psychosocial factors such as daily social activities and 
related psychological symptoms also provide a possible 
basis for the circadian variation in TTH. The temporal 
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relationship between headache intensity and psycho- 
social factors also needs to be investigated. It is of note 
that the distribution of headache exacerbations in TTH 
seems to be different than migraine. In most previous 
studies, migraine attacks were reported to occur more 
often in the morning [1,13,14] and sometimes during 
sleep [15], although a few studies showed peak attack 
frequency in the afternoon [16]. 

There are some limitations to the present study. 
Firstly, external clock time was used for the time axis in 
this study. This approach might cause between- 
individual variations in the effect of time, which might 
disappear with consideration of individual sleep-wake 
cycles. However, whether to use external clock time or 
individual sleep-wake cycles as the time axis is 
dependent on the hypothesis. Second, computerized 
EMA is not perfect, although it has an advantage over 
recalled self-reporting because it is able to capture the 
time course of headache intensity in a more detailed and 
accurate manner. Signal-contingent recordings were not 
continuous, but were sparse so that brief changes in 
headache intensity might have been missed, which was 
partly solved by dividing time into specific time blocks, 
sacrificing high time resolution. In addition, although 
correspondence between momentary headache intensity 
and recording time was accurate, there might have been 
a lag between the time of headache exacerbation and the 
time of making an event-contingent recording, despite 
instruction to make event-contingent recordings as soon 
as possible when headache was acutely exacerbated. 
Third, although we tried to explore the effect of medica- 
tion on the diurnal pattern of headache by analyzing the 
interaction with medication use, it would be better if we 
could investigate only patients who were free of medica- 
tion, especially for the investigation of pathophysiology. 
Fourth, the number of patients was rather small, thus 
the generalizability of this study might be limited. 

Conclusions 

Using computerized EMA, it was determined that there 
was a specific diurnal pattern of momentary headache 
intensity and the probability of acute headache exacerba- 
tions in TTH patients. These findings would be utilized 
for patient education and headache management. 
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